A multibaker map for shear flow and viscous heating 
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A consistent description of shear flow and the accompanied viscous heating as well the associated 
entropy balance is given in the framework of a deterministic dynamical system. A laminar shear flow 
is modeled by a Hamiltonian multibaker map which drives velocity and temperature fields. In an 
appropriate macroscopic limit one recovers the Navier-Stokes and heat conduction equations along 
with the associated entropy balance. This indicates that results of nonequilibrium thermodynamics 
can be described by means of an abstract, sufficiently chaotic and mixing dynamics. A thermostating 
algorithm can also be incorporated into this framework. 
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I. INTRODUCTION 

Shear flows provide one of the paradigms of transport 
processes The importance of chaos in the equa- 

tions underlying macroscopic shearing has recently been 
addressed by various numerical studies 1^-^, which to 
some extent were supported by kinetic theory and 
rigorous mathematical work |^jlo[|. In contrast, how- 
ever, a simple, exactly solvable model based on a low- 
dimensional chaotic dynamics — whose mixing property 
would be the cause of irreversibility — has not yet been 
established. For material and heat transport such type 
of models have helped to understand the physical con- 
tent of thermostating schemes used in numerical simu- 
lations (cf. however ^ for open questions). In 
the present article we introduce a similar model for shear 
flows in the hope that it can also serve such a purpose. 
The approach will be based on multibaker maps. Pre- 
vious works in this spirit successfully described the phe- 
nomena of diffusion [p^|-p^, conduction in an external 
fleld |p^|-^ , chemical reactions ||2^ , thermal conduction 
p^ , and cross effects due to the simultaneous presence 
of an external field and heat conduction 25 2^. 



Our aim is to model a sheared fluid confined between 
two parallel walls at the coordinates x — Q and x — L 
(Fig. |I]). The flow is assumed to be two-dimensional in 
the {x,y) plane. Shear is induced by prescribing differ- 
ent y components of the average velocities v of parti- 
cles close to the respective walls. In order to make the 
calculations more transparent, we conflne the discussion 
to cases where the driving is sufficiently weak to induce 
only a laminar flow, i.e., to cases where the velocity of 
particles is always directed in the vertical direction such 
that V = (0, ^(a;)). For this system we establish a local 
entropy balance that covers time dependent effects and 
does not rely on the implementation of boundary condi- 
tions. 

Three different boundary conditions for dealing with 
the dissipated heat are considered: (i) In the simplest 
case the system is isolated. A stationary linear velocity 



profile emerges in that case, and the temperature be- 
comes uniform. No steady state is reached, however, due 
to a constant increase of temperature in response to the 
viscous heating. In addition, we consider systems where 
(ii) there is a bulk thermostat uniformly taking out the 
viscous heat, and (iii) the temperature is fixed to the 
same value at both boundaries so that the asymptotic 
temperature profile is stationary, but no longer uniform. 
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FIG. 1. Graphical illustration of the shear flow. The sys- 
tem is conflned between two walls at positions x = and 
L. The walls move relative to each other in the y direction, 
thus inducing a velocity profile v{x) indicated by vertical ar- 
rows. For a flxed temperature at the walls this leads to an 
accompanying heat flow in the x direction (gray arrows). 

The paper is organized as follows: In Sect. || we recall 
basic notions of irreversible thermodynamics that are to 
be recovered in a suitable continuum limit of the multi- 
baker dynamics. In Sect. Ill we introduce the model, and 



establish the evolution equations for the velocity and the 
temperature fleld. This allows us to addres s th e entropy 
dynamics and its balance equation (Sect. |rv| ). Subse- 
quently, in Sect. ^ the macroscopic limit of the resulting 
equations is worked out. The global behavior at differ- 
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ent boundary conditi ons i s compared to irreversible ther- 
mod ynam ics in Sect. VII , and conclusions are drawn in 

Sect. Ivnl. 



II. IRREVERSIBLE THERMODYNAMICS 

In this section we recall the thermodynamic description 
of shear flows accompanied by viscous heating. The pic- 
ture is simplified by considering an incompressible fluid 
at constant pressure. 



A. Transport equations 

For a system with constant density and pressure the 
temperature T is the only relevant state variable, and for 
a complete description one also has to specify the veloc- 
ity field V of the fluid jl|,||. The thermodynamic state 
variables are the velocity field v and the temperature T 
1^ . Mass conservation is expressed by a continuity equa- 
tion. For incompressible fluids the uniform mass density 
p implies that the flow is divergence free, i.e.. 



diV,, = 0. 



(1) 



Here i = x,y labels the components of the local flow ve- 
locity V = (vxjVy), and we adopted the Einstein conven- 
tion, i.e., summation over repeated indices. The equa- 
tion of motion for the velocity components is given by 
the Navier-Stokes equation. For the case of a negligible 
pressure gradient it reads 



dt 



(2) 



where v is the kinematic viscosity, and d/dt is the total 
time derivative. 

The system of equations is closed [|| by the equation 

= —dAT -f ]- — {dkVi + divk){dkvi + divk) (3) 
dt pcv 2 cv 

for the temperature evolution. Here, cy is the specific 
heat at constant volume, and A is the thermal conduc- 
tivity. 



B. Entropy balance 

The balance equation of the entropy density s is 

dts = (t(*"' -I- $, (4) 

$ = (5) 

where the irreversible entropy production reflect- 

ing the viscous heating of the flow, $ the entropy flux, 
and j*^*^ is the entropy-current density. The term 



models an entropy flux let into a heat bath. In the bulk 
of a hydrodynamic system this term typically vanishes, 
but it takes non- vanishing values wherever there is a fiux 
into the environment (i.e., for instance at the boundaries 
when there is a heat flux through the walls). In simula- 
tions where one does not desire to focus on the effects of 
temperature gradients, and applies a thermostating al- 
gorithm, a non- vanishing can appear even in the 
bulk of the system. We say that a system is subjected to 
an ideal thermostat, when the dissipated heat is directly 
let into the heat bath. In that case the entropy current 
vanishes in the steady state, j^''^ = 0, while is non- 
zero in the bulk, and counterbalances the steady-state 
entropy production. 

For an incompressible fluid in local equilibrium the 
Gibbs relation ds = du/T applies locally. The evo- 
lution equations can be used to evaluate the terms in 
Eq. (Q). A straightforward calculation yields for the 
rate of irreversible entropy production (cf. for instance 
I, Ch. Xn. (23)]) 

(6) 

The associated entropy current takes the form 
Ch. XII. (22), (24)] 



-A 



T ' 



(7) 



It depends only on the local temperature and its gradi- 
ent. The flow velocity v does not enter. 



C. Laminar Flow 

For a laminar flow driven by prescribed non-trivial y 
components of the velocity at the two walls, the velocity 
field at any position {x,y) takes the form v = {Q,v{x)) 
(see Fig. |l|). The x component of the velocity vanishes, 
and the profile is translational invariant in the y direc- 
tion (parallel to the walls). We restrict our investigation 
to cases where the same holds for the temperature such 
that T — T{x). Consequently, for the laminar flow the 
transport equations take the form 



dtv — vd^v, 

d,T^—dlT+-{d.,v)\ 
pcv Cv 



(8a) 
(8b) 



while the rate of irreversible entropy production a^^^"^^ 
and the entropy current j^''-' can be written as 



(9a) 
(9b) 
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The scalar current j^^^ denotes the x component of the 
entropy current, and an analogous convention is adopted 
for all currents. The y components of the currents vanish 
in the considered setting. 
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FIG. 2. Graphical illustration of the action of the multi- 
baker map on its phase space {x,y). The indices labeling 
the columns and rows are given on the lower and left mar- 
gin, respectively, while their width is indicated at the up- 
per and right margins, (a) The mapping is defined on a 
domain of A'^ -|- 2 identical columns of square cells of size 
a X a that are labeled by the indices nix = 0, ■ ■ ■ , N + 1 and 
niy — — oo, • • • , oo. Boundary conditions on the fiow are im- 
plemented in the shaded cells tUx — and A^-l- 1. (b) The ac- 
tion of the map on regions that are mapped into cell {nix, nriy). 
The contraction and expansion for these regions is indicated 
by the deformation of the tags (R,U,S,D,L) in the different 
branches of the map. 



III. THE MULTIBAKER MAP 

In order to model the laminar flow we divide the plane 
into square cells of size a x a that are large enough, on 
the one hand, to admit thermodynamically meaningful 
averages, but, on the other hand, sufficiently small to ne- 
glect gradients across cells. The cells are considered as 
the regions used in irreversible thermodynamics to de- 
fine local equilibrium variables. The system shown in 
Fig. is represented by a rectangular array oi N = L/a 
cells in horizontal direction, and an infinite number of 



cells vertically (Fig. ga). Cells are labeled by the indices 
TTix = 0, • • • , iV -I- 1, and my = —oo, • • • , oo. All cells have 
the same dynamics except the outermost ones where it is 
modified to implement boundary conditions. 



A. Action of the mapping in the (a;, y) plane 

After each time unit r every cell is divided into five 
columns (Fig. ^). The rightmost (R) and leftmost (L) 
column have width ag. They are mapped onto a strip 
of height ag in cell (rUx + and cell (m^, — l,my), 

respectively. The left column (U) of widths au is mapped 
upwards onto a strip of height au in cell {rrix, my + 1), and 
the right one labeled by (D) downwards into (to^,, — 1), 
respectively. Region S stays within the cell. In all cases 
the area of the strips is preserved. 

This dynamics is driving two fields. Denoting the com- 
posite index {mx,my) as m, these fields are 

(i) the velocity field describing the mean flow velocity 
v{x,y) in cell m 

(ii) an energy field that represents the kinetic energy 
of the fluid in cell m. 

The dynamics on this two-dimensional lattice can be 
considered as a model of the velocity and energy trans- 
port in the conflguration space. In order to obtain a 
faithful representation of the entropy balance, however, 
one has to consider the phase-space dynamics. To define 
the dynamics in the analog of a /z-space, we take into 
account the translation invariance of the problem. The 
velocity and energy fields can take on different values in 
the columns rUx = 0, • • • , + 1 of Fig. ||a, but the fields 
have to be uniform within every column. In this respect 
it is not necessary to follow the dynamics in the y direc- 
tion, and one achieves a quasi one-dimensional dynam- 
ics in an {x,p) space, where p represents the non-trivial 
momentum-like variable of the model. 



B. Action of the mapping in the {x,p) space 

The domain of the multibaker in the {x,p) space is 
shown in Fig. It comprises a chain of -f 2 cells of 
size a X b, that for sake of more condensed notations are 
labeled by the index m. The middle N cells represent 
the bulk, and two additional ones are used to implement 
boundary conditions (cf. Section Vl). The parameter b 
sets the momentum scale. It will not play any role in 
thermodynamic considerations. 

In order to maintain the same dynamics in the trans- 
port direction in the two representations of the flow, each 
cell is divided into three columns of size ag, as and ag, 
where g is the same as above and s — l — 2g. The left and 
right columns of cell m are mapped into a strip of height 
ag in cell to + 1 and cell m — 1, respectively, as shown in 
Fig. 0b. The middle column of size as is squeezed and 
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stretched onto a strip of height as and remains in the 
same cell. 

This multibaker dynamics drives the velocity and en- 
ergy fields. Their values might depend on the phase- 
space coordinate. Hence, we are dealing with the bivari- 
ate distribution of the velocity v{x,p) and kinetic energy 
w{x,p) within each cell m. Only the cell averages, Vm and 
Wm appear in the transport equations. The dependence 
of ^(a;,^) and w{x,p) on the phase-space coordinates con- 
tributes essentially to the entropy dynamics, however. 
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FIG. 3. Graphical illustration of the {x,p) multibaker rep- 
resentation of the laminar flow, (a) Domain of the map. The 
outermost cells (indicated as grey boxes) are again used to 
implement boundary conditions, (b) Action in the bulk. The 
average values of the fields v{x,y) and w(x, y) on the ceUs 
[cf Eqs. and (^] and strips [cf Eqs. and (^] are 
given on the margins. 



Vm = {gp ) (12) 

at any instant of time. If the velocity is initially uniform 
in the full cell (i.e., the values of v in all strips agree), 
then due to Eqs. ( p^ ) and ( pT| ) the average velocity of cell 
m becomes after one time step 

P'mV'm = (1 ^ '^g)PmVm + 9Pm-lVm-l + gPm+lVm+l- 

(13) 

Observing that the mass density p is uniform and con- 
stant in time, the updated value of the cell velocity Vm 
is found to be 



(14) 



This evolution can be written in the form of a discrete 
balance equation 



Jm+l 



(v) 



with the discrete current 

,2 



(v) 



a g Vm - Vm-1 
T 



a 



(15) 



(16) 



In the boundary columns m — Q and iV + 1 the dynam- 
ics will be augmented by force terms (Sect. [v| ) in order 
to fix the velocity to the constant values = vl and 
UAT-f^i = tifl, respectively, irrespective of the momentum 
flowing into these cells. 



C. The dynamics of the velocity field 

The mass density pm is proportional to the average 
phase-space density p„j, in cell m, i.e., 

= Mhg.nr, (10) 

where M is a constant of dimension mass. When inter- 
preting the mass density or associated quantities from 
the point of view of the two-dimensional flow, p^ is to 
be understood as a mass density in the x coordinate per 
unit length in y direction. Hence, PmVm represents the 
preserved momentum density (per unit length) of the 
hydrodynamic flow. Its evolution expresses momentum 
conservation. As a consequence, the updated values for 
the velocity on the strips (L, 5, R) after one time unit are 
(cf. the right side of Fig. yb) 

P'mV'm.l = Pm+lVm+1, 

P'mV'm.s = PmVm, (H) 
PmVm.r ~ Pm — lVm — lj 

where, v'^ v'^ ,., and v'^ ^, stand for the velocities of the 

flow in the regions (L, S, R) of cell m, respectively. The 
prime indicates that the updated values of the fleld are 
considered. The average momentum PmVm of the full cell 
is the average of the contributions on the different strips, 
i.e., 



D. The dynamics of the energy fleld 

The energy of cell TTi is obtained by integrating CmQm 
over the volume ab of the cell. At any instant of time it 
is the spatial average of the energies 6^,1, &m,s and em,r 
of the different strips. 

The difference between the energy e and the translational 
specific kinetic energy Mv'^/2 of the flow deflnes the spe- 
cific kinetic-energy density w = e — whose macro- 
scopic limit will be proportional to the local temperature. 
Therefore, on the strip i = l^s,r in cell m one observes 
kinetic-energy densities Wm,i that fulfill 

_ M 2 

Gm,i — ~^Vjy^ ^ -\- Wm,i, 

while the coarse-grained kinetic energy Cm obeys at the 
same time 

M , 

em = Y^m+V'm- 

Using these definitions and the averaging rule ( p^ ) for 
the energy one finds 



M 

Wm = 9Wm,l + SWm,s + gWm,r + 



x2 



(18) 



This shows that the average kinetic energy Wm is not the 
straightforward spatial average of the values Wm.i on the 
strips. Rather intra-cell variations of the velocity field 
(i.e., a non-trivial distribution of the Um.i) also contribute 

to Wjn- 

For a thermally closed system the specific kinetic en- 
ergy w is advected by the flow. The values on the strips 
after one time unit are then 



(19) 



From this and the update ( |11| ) of the velocities the up- 
dated kinetic energy is found to be 

w!^^Wm+ .g(Wm-l + W„i+i - 2w„) -|- 



+ 



2 

M ' 

T 



{Vm-l 



+ {Vm+1 - VjnY 
2 



a 9 



(20) 



The terms proportional to Af describe the effect of vis- 
cous heating in this discrete setting. 

In order to model the action of a thermostat, that leads 
to local changes of the specific kinetic energy due to a 
heat flux into the environment, an additional thermostat 
heat source q-m is incorporated into the update of Wm by 
multiplying the right-hand side of Eq. ( po| ) by a factor 

[1 + Tgm], 



+ 9{w,n^i + W^+i - 2w„i) + 



gM 

2 
M 



{Vrn-l - Vmf + [Vm+l " ^m)^ 



9 Ki-i + v„i+i - 2w„ 



This equation can be rewritten as 

W'^-Wm _ ^ ^ a^9 Wm+l - 2w„i + W.m-1 
T T 

with the full heat source 



(21) 



(22) 



1 + Tq„ 
a'^9 M 



'Jm+1 



2 



a" 

-I 2 



(23) 



The equation for the update of w can be rearranged into 
a balance equation for the heat per unit volume 



.{w) _ .{w) 
Jm+1 Jin- 



(24) 



which comprises the divergence of the discrete "heat" 
current 



^"^ ~ T M a 



(25) 



The first contribution to Qm in ( |23| ) refiects the action of 
the thermostat, and the latter two the effect of viscous 
heating of the fluid. A steady state with a uniform w pro- 
file can be found for Q ~ 0, which thus mimics an ideal 
thermostat. In contrast, in a bulk system in the sense 
of conventional irreversible thermodynamics q^ ^ 0, and 
Qm only vanishes when the discrete velocity gradients 
Vm+i — Vm and Vm — Vm-l across the cell vanish. 



IV. ENTROPIES AND THEIR TIME 
EVOLUTION 

A. The coarse-grained and the Gibbs entropy 

The Gibbs entropy for cell to of a multibaker system 
with a density field q{x,p) and specific kinetic-energy 
field w{x,p) takes the form p5|] 



sin ^ 



cell 



dx dp q{x,p) In 



q{x,p) 



w{x,py 



(26) 



where 7 is a free exponent. Since the phase-space density 
q(x,p) is constant, one can identify the reference density 
Q* with q{x,p)^ and obtains 



M 



cell 



^fln^x,p). (27) 
a 



Here, we have already replaced the phase-space density 
by the mass density (p^. 

The coarse-grained entropy of cell m is based on the 
cell-averaged coarse-grained value Wm of the kinetic en- 
ergy w{x,p). For the form ( p7| ) of the Gibbs entropy, one 
obtains 



^sa-f— In 



(28) 



as the coarse-grained entropy. 

In order to discuss the thermodynamic time evolution 
of entropies, one conveniently starts with uniform densi- 
ties in every cell [^,|l^. In that case the coarse-grained 
and the Gibbs entropy initially coincide. After one time 
step, however, they typically differ. The Gibbs entropy 
has changed due to the fact that the w field takes dif- 
ferent values on the strips {R,S,L), leading to the new 
value 
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S^^^ = kB-fgab {g ln[it;„_i (1 + rq,„)] + (1 - 2.g) ln[w,„ (1 + rq,„)] + g ln[it;,„_i (1 + rq,„)]} 
= ks'-fgab \ \n[wm (1 + rq™)] + 5 In - g In ■ ^" 



Wm+l 



(29) 



On the other hand, the coarse-grained entropy after one 
time step is 



5^ — —kBCib'-fglnw'^ 



(30) 



and 



T M Wm 



(35a) 



B. Entropy balance 



(35b) 



The coarse-grained entropy only depends on the aver- 
age kinetic-energy density Wm, and is therefore consid- 
ered as the multibaker analogue of the thermodynamic 
entropy. Its temporal change can be decomposed as 



describes the direct flux into the thermostat. Such a flux 
is encountered whenever there is a non-vanishing source 
term qm- In view of Eq. ( p3| ) this finding further sup- 
ports th e inte rpretation of q,n and Qm given at the end 
of Sect. FTTdI. 



abr abr 

(^Sm Srn ^ ) (^Srra Srn ^) 

abr 



a{G)' _ a(G) 
abr 



, (31) 



V. THE MACROSCOPIC LIMIT 



where information-theoretic arguments [ pO| , pl 25 1^ lead 
one to identify 



abr abr 
with the rate of entropy production, and 

oiG)' _ oiG) 



abr 



abr 



(32a) 



(32b) 



with the entropy flux. Note that the second term of the 
numerator of /S.iSm vanishes due to the initial condition 
of uniform fields in the cells. 

In serting the values for S'^ and Sm into the definition 
(HI) yields 



ar 



kslP 



Mt 

Wm-l 



ln^(l + rg„0^i 

Wm 



(33) 



gin 



gin 



Remarkably, this expression for the entropy production 
does not depend on the source term Qm, but only on 
the values of the coarse-grained field w in cell m and its 
neighbors. The shear flow enters only indirectly through 
the update w'^^ of Wm- 

Similarly to the other balance equations, the entropy 
flux can be written as 



ar 



Jm+l Jm 



(34) 



(s) 

where the discrete entropy current jm takes the form 



In this section we evaluate the expressions of the differ- 
ent quantities considered in Sects. Ill and [V, and iden- 
tify conditions for consistency with the thermodynamic 
results described in Sect. ||. 



A. Definition of the limit 



The macroscopic limit implies that L ^ a (i.e., N 3> 
1), and T is much smaller than typical macroscopic time 
scales. Formally it can be deflned as 



a, T 







such that the spatial coordinate 



(36) 



(37) 



is flnite. As mentioned earlier, the fleld w is assumed to 
go over into the local temperature T[x) in the macro- 
scopic limit, i.e., 



Wm CkBT{x), 
where C is a dimcnsionlcss constant. 



(38) 
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B. The transport equations 

The macroscopic form of the velocity cmrent (|l^) be- 



comes 



9 



(39) 



In order to achieve a meaningful thermodynamic result 
the ratio a^g/r has to be finite in the limit. Indeed, 
comparison with Eq. (|8^) shows that 



a^g 



(40) 



is the kinematic viscosity. 

Equations (|3|) and (|^) imply that Eq. (O) is the 
discrete form of the Navier-Stokes equation ( |8a| ) for the 
laminar flow. Moreover, with this choice for g one obtains 
for the heat current (P5| ) 



(41) 



Hence, Pick's law of heat conduction is recovered with 
the heat conduction coefficient 



(42) 



Due to Eq. 
duces to Eq. I 

dtT^ 



3) , the macroscopic limit of Eq. ( |2l| ) re- 
^ for the temperature evolution: 



A 



CkBp/M 



C/cr 



qT. 



(43) 



Comparing the coefficients in this equation with the ones 
in Eq. (pq), one obtains 



CkB = cvM. 



(44) 



The proportionality constant C introduced in Eq. ( |3q ) 
corresponds thus to the specific heat at constant volume 
(measured in units of fcs). 



C. The entropy balance 



In the macroscopic limit the rate of irreversible entropy 



production (33) becomes 



(irr) 



(45) 



It fully agrees with the thermodynamic form of the en- 
tropy production ( pa] ) when the coefficient k^jvp/M in 
front of the first term is the heat conductivity. In view 
of (42) we thus conclude, that in our model 7 = C, i.e., 
the exponent 7 appearing in the definition (|2^) of the 
entropy is proportional to the specific heat. Hence, the 
final form of the heat conductivity A can be settled to 



pjkBV 
X = ^^^pcvv. 



(46) 



Working out the expression (35a) of the entropy cur- 
rent jm\ one finds in the macroscopic limit 

j(^){x) = ~kBiy—— = -\—. (47) 

In view of (^) this relation also fully agrees with its 
thermodynamic counterpart (9b). 



VI. BOUNDARY CONDITIONS 

In this section we demonstrate how boundary condi- 
tions for the shear flow can be implemented in the multi- 
baker dynamics. 
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FIG. 4. Implementation of the boundary condition in cell 
0. The indication of the width of strips, values of the fields 
on the strips, and action of the map follows the rules spelled 
out in Fig. ^ 



A. Boundary conditions for the velocity field 

The multibaker dynamics for cell is shown in Fig. IJ. 
(An analogous prescription holds in cell + 1.) The ac- 
tion of the strip mapped from cell to 1 fully agrees with 
the one of the corresponding strips in the bulk such that 
the dynamics of cell 1 agrees with the bulk dynamics of 
Fig. gb. The important difference about the dynamics at 
the boundary is that 

(i) there is no column leaving cell to the left, and 

(ii) the momentum of the particles entering cell is not 
preserved. 

Condition (i) reflects that no particles can penetrate the 
walls of the system, and (ii) reflects that there is a force 
exerted by the boundary on the fluid. In strip L entering 
cell (and analogously strip R entering cell A^ -I- 1) the 
update of the velocity consequently contains the contri- 
bution of the external force /o (/at+i), leading to 



V 



N+l,r 



= + 



9P ' 
Jn+it 

9P 



(48) 
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The update of the coarse-grained velocity of the leftmost 
cell is 



Q 



(49) 



and an analogous relation holds for cell + 1 . A shear 
flow with a prescribed shear rate is enforced by adjust- 
ing the forces in such a way that uq and wat+i take the 
respective values Vh and vr of the velocities at the walls 
at any time. 

In view of the time evolution of Wm, in a steady state 
[Eq. (|lj) with v'^ = Vm\, the asymptotic velocity profile 
becomes linear 



Tfl TTL 

{vr - vl) = wl + -TTTT (50) 



N +1 



N + 1 



irrespective of the temperature distribution in the sys- 
tem. 



B. Boundary conditions for the kinetic-energy 
density 

The modification of the dynamics of the velocities at 
the boundary also implies changes of the thermostat 
source terms. Since the velocity in cell is uniform by 
definition, the update of wq can be found from Eq. jlq ) 
where the terms containing Wm,i and Vm.i are not present 
since no particles enter cell from the left. Consequently, 

w'q = [wo + g{wi - Wo)] (1 + Tqo). 

This leads to the discrete temporal change 

^2 



w, 







Wo 



90 



1 



"0 



1 a g wi 
a T 



7«0 



(51) 



instead of Eq. (p2[). For a steady state w'q — wq vanishes, 
and the equation implies [cf. Eq. (kF 



ago 



{wi - Wo)/a 
wq + g{wi - Wo) 



Ji 



g [wq + g{wi - Wo)] ' 
(52) 



The right-hand side of this expression is finite in the 
macroscopic limit, where 



ago 



T 



AT(0) ■ 



(53) 



Since, however, a ^ in the macroscopic limit, the ther- 
mostat heat source qo cannot be interpreted as a density 
when there is a finite heat current j^'"-'(0). Its integral 
over the cell, aqo, stays however finite. As a consequence, 
the heat current has to vanish in a steady state when 
(Jo = qN+i = 0. 



VII. DIFFERENT MACROSCOPIC FLOWS 

We have seen that in the macroscopic limit the local 
momentum, energy and entropy balances for the multi- 
baker map coincide with their thermodynamic forms. 
This result was achieved by inspecting the local time 
evolution of the densities without referring to particular 
boundary conditions, and it is not restricted to steady 
states. In order to underline these features, we discuss 
the global transport pictures for three different settings 
of transport. We emphasize, that the considerations im- 
mediately generalize to arbitrary time-dependent states 
although the particular calculations are carried out for 
stationary coarse-grained velocity fields. 



A. Isolated systems 

In thermodynamics sheared systems are often consid- 
ered to be thermally isolated. In that case q,„ identically 
vanishes for all cells to = 0, • ■ • , A'' -I- 1. Eq. ( |5l| ) implies 
for the current (Eq) 



Ji — 



a^g p wi- Wo 



ap Wq — Wo 



(54) 
fi- 



T M a M T 

Assuming that the time derivative of wq is 
nite in the macroscopic limit, the current jj""' 
—a{p/M)CkBdtT{0) approaches zero as a ^ 0. An 

analogous arguments applies also to j^^. 

Thus, for the asymptotic state, the boundary condi- 
tions imply that wo and wn+i practically coincide with 
wi and wn, respectively. Moreover, the viscous heat- 
ing term Qm in the bulk only depends on the square of 
Vm+i — Vm- The heating is therefore spatially uniform 
in the long time limit. With this input, one immediately 
verifies that w,„ approaches a spatially uniform value w* . 
However, in view of the heat source (5,„ in (p^), it con- 
stantly grows in time, 

w* — w* Av \^ 

= i/Af — . (55) 

This temporal evolution of w reflects the rise of tempera- 
ture due to the entropy production in the system. After 
all. 



In- 



kBlP 



Av 



0{r) 



tAI w* w* ' \a{N + 1) 

which in the macroscopic limit reduces to the thermody- 
namically expected value 



T*{t) 



Av 

T 



where T*{t) — const -|- {v / cv){Av / L) t is the spatially 
uniform, but temporally increasing asymptotic tempera- 
ture distribution. Note that the flux and the en- 
tropy current vanishes in this setting, implying an ever 
increasing entropy. 
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B. Systems subjected to an ideal (bulk) thermostat 

In a system subjected to an ideal thermostat the vis- 
cous heating is instantaneously released into a heat bath. 
Accordingly, Qm vanishes for all cells m — 1, - ■ ■ ,N such 
that Qm can be determined from Eq. (|2^). 

M oreove r, according to the results obtained in 
Sect. VI B we are again dealing with a system with van- 
ishing heat currents through its boundaries, i.e., for the 
asymptotic state the boundary conditions imply wq = wi 
and wn = wn+i, and there is no source term Q for the 
heat. Consequently, Wm approaches a spatially uniform 
state w* that is stationary in this case. Entropy produc- 
tion arises due to the nontrivial form of qm [cf. Eq. (|2^ 
with Qm = 0] 



( At) 
Vq(W+1) 



vMt 



( At) 

\a{N+l] 



in the bulk [cf. the steady-state velocity profile 
Hence, the irreversible entropy change is 



ln(l + Tg*) 



V 



Mw* a^{N + iy 



Oir), 



and in the macroscopic limit it reduces to the thermody- 
namically expected value 



releasing the viscous heat into the bath at the bound- 
aries [cf. Eq. (|5^)]. In contrast to the previous cases, the 
discrete heat current does not vanish at the boundaries. 
In view of Eqs. (|2^) and (|5l| ) the heat current through 
the left boundary of the steady state system is 



] 



(w) 



p 



Wi - Wq 



a 



ago 



1 + rqo 



Wq. 



(58) 



Since the dynamics is symmetrical and wn+i = wq, the 
current at the right boundary takes the same value up to 
a change of the sign 



.(«.) 



Jl ■ 



In a steady state the cells in the bulk consequently fulfill 



= p- 



Jm+l 



r ■ a2(7V + l)2 

leading to a linear profile of the heat current. 



(f) = a I TO - 1 



N 



pv 



Av 



a{N+l) 



In view of the definition Eq. (^5|) of jm \ this implies a 



(59) 



(56) 

parabolic profile of w in the steady state 



Wq a m - (N + 1) M — 



, (60) 



(irr) 



pv f All 



T* \ L 



where T* is the steady-state temperature. Also in this 
case there is no entropy current in the steady state. How- 
ever, at every location in the system there is a non- 
vanishing entropy flux into the thermostat 



Av 
L 



(57) 



which exactly compensates cr^*'''") such that the entropy 
indeed is stationarity. 



C. Thermostating at the walls 

Finally, we discuss a steady state in a hydrodynamic 
bulk system that generates heat flux into the walls due 
to the prescribed temperatures at the boundaries. There 
are no thermostat heat sources in the bulk, i.e., q„i = 
for TO = 1, • ■ • , iV. On the other hand, there are source 
terms go and qN+i — qo acting in the two outermost 
cells that fix the values of wq and wn+i (i.e., the tem- 
perature) to the same constant value wq- To this end 
the sources counterbalance a macroscopic heat current 



as expected for the temperature profile in a shear flow 
subjected to identical temperatures at the two ends. 

The entropy production in the bulk is related to the 
spatial variation of the steady-state distribution w* 
[cf. Eq. (132 



(irr) — 



' tM 



gin 



w„ 



w 



gin 



M 



1 _ of"") 

w* a 



Jm+l 



' Jm 



2pvw. 



*2 



, (61) 



where the second equation was obtained by expand- 
ing the logarithms to second order in the differences 
-iwl+i — wt, , and rearranging terms. Using Eqs. (|24|), 
(P3| ) and ( |59| ) , the first term can be related to the applied 
shear rate, while the second one represents the contribu- 
tion to the entropy production arising from the heat flow. 
Indeed, in the macroscopic limit one recovers in this case 
both contributions to Eq. (^5|). (For small shear rates 
the contribution from the temperature change is, how- 
ever, negligible). Correspondingly, there is a finite en- 
tropy current j*-"' at every point in the system, but only 
at the boundaries there is a flux into the thermostat. 
According to Eq. (^^ the macroscopic limit of the full 
entropy flux through the boundaries becomes 
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M T 



(62) 



Since the stationary temperature profile is obtained 
from @) in the form of T{x) = T(0) - [x {x - 
L){M/2CkB)]iAv/Lf, the derivative at the origin IS 
9xT\^^Q = ML{Av/L)'^/{2CkB). The contribution at 
the right boundary is the same. Therefore, we find for 
the integrated entropy flux 



'To' 



L 



Av 

T 



(63) 



where we used 7 = C. For weak shear this is essentially 
the same as the integral of the constant flux ( |57| ) over the 
chain. Thus, we conclude that for sufficiently small heat 
currents thermostating in the bulk and in the boundaries 
can lead to the same global behavior. 

The work per unit time done by the external force 
densities /o and Jn+i can be evaluated using their form 
Eq. (B3) and the value of the constant velocity gradient 



aifoVL + In+iVr) = vp 



[Avf 



(64) 



This expression is exactly — Tq^'^*'''*"*-', such that in a 
steady state the work done by the external forces equals 
the total heat flux into the thermostat. 



VIII. DISCUSSION 

We have enlarged the family of multibaker maps by a 
model for momentum, energy and entropy transport in 
viscous hydrodynamic flows. Although the macroscopic 
problem is of strongly dissipative nature, the proposed 
multibaker dynamics is area preserving (Hamiltonian) . 
This is to be contrasted with previous multibaker models 
of electric transport |^,^ and thermoelectric cross ef- 
fects 1^,^ where the inclusion of a reversible dissipation 
mechanism was necessary in order to simulate the effect 
of thermostating on the particle dynamics. The form in 
which thermostating appears in the present model is via 
a heat-source term in the microscopic energy dynamics 
(which was already present in the model of cross effects 
p5|). This term, however, does not give rise to phase- 
space contraction. 

The model has the following basic features: 

(a) The time evolution of the system can be interpreted 
as that of weakly interacting particles. The resulting 
"multi-baker" gas obeys the classical ideal-gas equation 
of state. The particle and heat diffusion coefficient, and 
the kinematic viscosity are proportional to each other, as 
in the kinetic theory for classical ideal gases . 

(b) The distribution of a macroscopic velocity does not 
enter the entropy explicitly. Rather the shear rate ap- 
pears in the entropy balance via its influence on the tem- 
perature dynamics only. 



(c) The connection to a thermodynamic description of 
transport is achieved by considering fields, which are 
coarse grained in regions of small spatial extension. Their 
properties are to be contrasted to those of fields charac- 
terizing the microscopic evolution. 

(d) Comparing these two levels of description allows us 
to identify all contributions to the local entropy balance, 
in full consistency with thermodynamics. They apply to 
both stationary and transient states. 

(e) A source term is implemented in the evolution equa- 
tion of the microscopic kinetic energy. It provides the 
possibility to implement local irreversible cooling of the 
system, i.e., to extract heat such that states which are 
permanently warmed up by viscous heating can become 
stationary. In traditional thermodynamics these terms 
are only present at the boundaries, and they vanish in 
the bulk. 

(f) It is indicated how the velocity and the kinetic- 
energy dynamics can be implemented at the walls in or- 
der to achieve correspondence with different macroscopic 
boundary conditions. 

(g) When source terms are present in the bulk, the lo- 
cal entropy balance of nonequilibrium thermodynamics is 
generalized by introducing at every location an instanta- 
neous flux of entropy (i.e., of heat) into a thermostat (the 
entropy flux is then no longer the divergence of the en- 
tropy current). The dynamics is in that case reminiscent 
of numerical algorithms based on Gaussian thermostats. 

(h) The global entropy balance was worked out, in order 
to demonstrate that the total heat flux into a thermostat 
is independent to a large extend of whether thermostat- 
ing is applied in the bulk or at the walls. 

The major interest of the present model lies in the light 
it sheds on the origin of viscous heating in determinis- 
tic models of transport. It was pointed out how fractal 
structures emerge in multibaker models for a variety of 
physical settings of shear flow. In all cases the struc- 
tures arise from the mixing of regions with different local 
temperatures and flow velocities whose differences are ex- 
ponentially proliferating to smaller and smaller scales for 
a nonequilibrium system. One can explicitly follow this 
redistribution of the kinetic energy, until it reaches the 
scale of the domains used to define local thermodynamic 
averages. Motion at that scale has to be considered as 
contributing to the non-directional motion, hence leading 
to viscous heating. It is only due to this coarse-graining 
mechanism that the macroscopic shear rate appears in 
the expression of the irreversible entropy production. 
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